Coherent Dynamics of Vortex Formation in Trapped Bose-Einstein Condensates 
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Simulations of a rotationally stirred condensate show that a regime of simple behaviour occurs in 
which a single vortex cycles in and out of the condensate. We present a simple quantitative model 
of this behaviour, which accurately describes the full vortex dynamics, including a critical angular 
speed of stirring for vortex formation. A method for experimentally preparing a condensate in a 
central vortex state is suggested. 
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The production of vortices has been a central issue in 
the study of superfluids. It has been demonstrated for ex- 
ample that attempts to produce bulk rotation in a cylin- 
der of He II will lead to vortex production, a state which 
calculations show to be energetically favoured. The cur- 
rently realised ||l|,|[ Bose-Einstein condensates (BEC) of- 
fer a new medium for studying vortices, and a number 
of theoretical studies have considered the properties of 
static vortices |3[^, their stability ||-|§1, excitation spec- 
tra ||,|l^, and phase sensitive detection techniques pi[ . 
A variety of methods have been considered by which vor- 
tices could be formed in a BEC. Fetter g has suggested 
that a rotating nonaxially symmetric trap could imitate 
the He II rotating cylinder, and obtained an approxi- 
mate critical rotational speed for vortex production by a 
heuristic argument. Jackson et al. [n2| have shown that 
vortices may be generated by movement of a localised po- 
tential through a condensate, while Marzlin and Zhang 
p3[ have calculated vortex production using four laser 
beams in a ring configuration. Other numerical simula- 
tions, such as collisions of condensates H,|lHl have shown 
in fact that vortex production appears to be a common 
consequence of mechanically disturbing a condensate. 

In this paper, we consider a trapped BEC stirred ro- 
tationally by an external potential, and find and analyse 
a regime where only a single vortex forms. We present 
a simple quantitative model of this behaviour, which ac- 
curately describes the full vortex dynamics in terms of a 
coherent process. The model gives the critical speed of 
rotation for vortex formation and explains a number of 
other features that are seen, including the stability of a 
central vortex (at T = 0). Our investigation is based on 
the Gross-Pitaevskii (GP) equation for the condensate 
wavef unction ^{r,t) 



. di>{v,t) 

' dt 



+ C|V(r,t)|V(r,t), 



(1) 



which is known to accurately describe condensates close 
to T = 0. In Eq. (|l|) we have used scaling and notation 
as in Ruprecht et al. Ilffl, V{r,t) is the external poten- 



tial, and C is proportional to the number of atoms in 
the condensate and the scattering length. We consider 
the GP equation in two dimensions only, and solve it nu- 
merically. We simulate the effect of stirring by adding to 
the stationary trap potential a narrow, moving Gaussian 
potential, representing, for example, a far-blue-detuned 
laser Q. V{r,t) is given by r^/4 -|- W{r,t), and the 
stirring potential 



W{r,t) = H/oexp 



|r-r.(i)| 



Ws 



(2) 



is centred at rs{t). The initial condensate state for our 
simulations is the lowest energy eigenstate of the time- 
independent GP equation M in which V includes the 
stationary stirrer. The stirrer moves anticlockwise on a 
circular path, accelerating constantly until i = tt, when 
it reaches its final angular speed uif. Figure ^, which 
shows the state of a condensate after it has been stirred 
for some time {t < Stt) then left to freely evolve {t > Stt) , 
illustrates the complexity of behaviour that can occur. 
Vortices of positive and negative circulation have formed 
and, as time progresses, move relative to each other, and 
annihilate when a positive and negative pair collide [[l8[ . 
Amidst the complexity of possible behaviours, an im- 
portant and simply characterised behaviour emerges, 
namely the formation and dynamics of a single vortex. 
An example which illustrates the main features is given 
in Fig. 2 where sequential subfigures show the evolution 
of the condensate as the stirrer revolves. A single vortex 
enters at the edge of the visible region of the condensate, 
then cycles to the centre of the condensate, and back to 
the edge. This cycle repeats regularly, as can be seen in 
Fig. where the solid line shows the angular momentum 
(L) plotted as a function of time for this case. At lower 
stirring speeds, similar vortex cycling occurs, but with 
progressively smaller amplitudes as cuf decreases, so that 
the vortex oscillates near the condensate edge. We have 
found that the condensate gains angular momentum even 
for very small values of ujf. The critical angular speed 
which causes a single vortex to cycle right to the centre 
of the condensate we denote ujc- In Table | we present 




FIG. 1. Probability density at t = 127r of a condensate 
stirred as described in the text, with the stirrer gradually 
withdrawn between t = Att and 5tt. Contours are logarithmi- 
cally spaced. Vortices are detected numerically by searching 
for their 27r phase signature, and are marked near dense re- 
gions of the condensate by a + or — sign according to their 
sense. C = 88.13, ujf = 1, r, = 3, Wo = 10, Ws = 1. 

results from our simulations that show that lUc decreases 
as C increases, in agreement with the heuristic result of 
Fetter §. 

The single- vortex behaviour can be understood quanti- 
tatively in terms of a nonlinear Rabi cycling model. The 
essential idea is that the stirring potential causes the con- 
densate to cycle between the ground state and the first 
vortex state, analogous to the Rabi cycling of an atom 
in a light field. We decompose the condensate on a ba- 
sis of a ground-state-like part (axially symmetric) and 
a vortex part (axially symmetric with an anticlockwise 
phase circulation). In the linear (i.e. C — 0) limit the 
condensate would be represented as a superposition of 
the ground state and the first vortex state of the trap. In 
the nonlinear system, it is more accurate to decompose 
the system into coUisionally coupled states, in which the 
radial form of each of the basis states is modified by its 
collisional interaction with the other. Accordingly we as- 
sume that the condensate mean-field wavefunction can 
be represented approximately as 



tp{r,t) = as{t)4>s{r,ny) + ay(t)(j)y{r,ny)e^ 



(3) 



where r and 9 are the cylindrical polar components of 
r, and Uy = |ai,p. We obtain the lowest energy coupled 
eigenstates (j)s{r,ny) and </)„(r, n„)e*^ together with their 
eigenvalues Hsiny) and ^y{ny) by solving for a particu- 
lar value of ny the coupled time-independent radial GP 
equations. 
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FIG. 2. Sequence of states for a condensate stirred from 
rest as described in the text. Probability density is shown in 
the left-hand column, and the phase of ^j) in the right hand 
column, for (a) t = 8.80, (b) t = 18.35, and (c) t = 28.15. The 
circle denotes the stirrer. Parameters are as in Fig. hi except 
ujf = 0.5 and Ts = 1.5. 



+ C {n„4>l + 2nA0^) 



(4) 



Here a and A are either s or u, /^ = and ly = \ 
are the angular momenta of (ps and 01, e*^ respectively, 
Uy is the fraction of the condensate in the vortex com- 
ponent, and ris = 1 — riy the fraction in the sym- 
metric component. The (j)„ are real nonnegative func- 
tions normalized as J (j)^dr — 1, and (j)s and (pye''^ are 
of course orthogonal. The superposition in Eq. (pi) pro- 
duces a condensate with angular momentum expecta- 
tion value (L) ~ Uy , and a vortex whose distance 
from the centre of the trap decreases as n„ ^ 1. In 
the absence of a stirrer, the vortex precesses about 
the centre of the condensate at a frequency [ly — ^s- 
Substituting Eq. (p|) into Eq. (nh, and projecting al- 
ternately onto the states (f)s and (pye''^, we obtain a 
pair of coupled equations for dg and dy. Noting that 
a constantly rotating stirring potential W{r,t) can be 
written e-''^^*^M^'(r)e+*"^*^, and writing a^ = a^e*"', 
dy = aye''^"^+'^f*\ where as(t) = J^^ fisit')dt' , we collect 
the oscillating exponential time dependences and trans- 
form to a frame which rotates with the stirring potential 
to obtain the equations 



— -f. = -iSsiny)ds - -fl{ny)dy, 



(5a) 



(L>1 




FIG. 3. Angular momentum expectation values versus time 
for ujf = 0.5 (solid line), ujf = 0.4 (lower dashed line), 
ujf — 0.6 (upper dashed line). Other parameters as in Fig. 2. 



—J- = -i[A{ny) + 6y{ny)]dy - -fl*{ny)ds 
Here A(n„) — /it,(n„) — /is(n„) — ujf and, 
Scriny) = / (t)ainy)W' {r)(t)cr{ny)dr, 
n{ny) = 2 j dps{ny)W'{v)(j)y{ny)e'Ur. 



(5b) 

(6a) 
(6b) 



Eqs. (0) formally resemble the classic Rabi equations, 
and hence we identify the J^ as frequency shifts and O as 
the bare Rabi frequency, but note that here these quanti- 
ties are variable and depend on the value of n„. Despite 
this nonlinear dependence, the concept of Rabi cycling 
provides a simple framework in which to understand the 
formation and dynamics of a single vortex: the stirring 
potential couples and causes cycling between the initial 
ground state and the first excited vortex state. The en- 
ergy E' in the frame rotating with the stirrer (obtained 
from the expectation value oi H' = H — UfL where H 
is the lab frame Hamiltonian and L is the dimensionless 
angular momentum operator) is conserved, and thus in 
any solution to Eqs. dq), hg and hy must follow a trajec- 
tory that conserves E' . Complete cycling of the vortex 
to the centre of the condensate occurs when Uy reaches 
the value of 1, but this requires the energies in the rotat- 
ing frame of the ground state (f>s{r,ny = 0) and the first 
excited vortex state (l)y{r,ny — l)e'^ to be equal. Thus 
the critical angular speed Wc is given by the relation 



Ey 



E. 



9' 



(7) 



where Eg and Ey are the lab frame energies of the ground 
state and first excited vortex state respectively. A finite 
stirrer shifts these energies by 5s{ny = 0) and 5y{ny — 1) 
respectively, adjusting Wc by their difference. 
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6.079 
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18.860 
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TABLE I. Critical angular frequency lOc for the 
two-dimensional condensate. The final column gives bounds 
for UJc found from our simulations of the full GP equation. 



In Table | we list the critical angular speeds predicted 
by Eq. (0) for a range of C values, along with the values 
of LOc found from our numerical simulations of the full GP 
equation for C — ZQ and C — 88.13 cases. The agreement 
between the predictions from the two-state model and the 
full numerical simulation is excellent. 

It is difficult to obtain accurate simulations at large 
values of C for numerical reasons. The C = case is 
easily tractable, and although no visible vortex cycling 
occurs below LOf — 0.25, a multiple vortex regime is en- 
tered at LOf = 0.8 < uJc\ the reason being, as Marzlin and 
Zhang |1^ have noted, that the trap levels are equally 
spaced for the linear case, so that mixing to higher vor- 
tex states readily occurs, and our two state model is no 
longer valid. The success of the two state model is de- 
pendent on the fact that for C 7^ the spacing of the 
levels is nonuniform. 

The Rabi model also allows us to explain other features 
of the behaviour, such as the period of cycling, the re- 
sponse to smaller stirring speeds, and the effect of differ- 
ent values of stirring radius rg ■ In Fig. ra the condensate 
response to stirring just below the critical speed is shown, 
and reveals an increase in oscillation frequency and de- 
creased transfer to the pure vortex state, compared to 
the critical case. In a simple Rabi model, where the de- 
tuning A and Rabi frequency fi are constant, the cycling 
frequency is fi = vfP + A^, and the maximum value of 
Uy is [Vt/Vl Y . By identifying the effective detuning for 
the two-state system to be A -I- 15„ — (5g, these expressions 
give a qualitative description of the subcritical stirring 
in Fig. 0. A more quantitative treatment requires the 
nonlinear character of Eqs. (0) to be taken into account, 
which is achieved by solving the coupled pair in Eqs. (0) 
to find the eigenvectors and eigenvalues at each value of 
Uy and then using these to solve Eqs. (0). We note that 
in Eqs. (0) the term 2n\(j)'j^ gives rise to an energy bar- 
rier between the Uy — and Uy = 1 states of the system. 
The constraint on the system imposed by the ansatz of 
Eq. (ph increases this energy barrier slightly compared to 
the true (unconstrained) case, and the accuracy of our 
procedure can be improved by decreasing this factor of 
2. For example, at C = 88.13, if we decrease the factor of 
2 — > 1.58, our two-state model produces behaviour which 
closely matches the results from the full GP equation, as 
we show in Fig. El The energy barrier is deformed by 
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FIG. 4. Angular momentuin versus time for the full GP 
equation simulation (solid line) and the two-state model 
(dashed line). The two-state model starts at t = 3 with 
ha = \/0.965 and a„ — — \/0.035. Parameters as in Fig. 2. 



the presence of the stirrer, allowing the system to cycle 
between the vortex and ground state. If the stirrer is far 
from the centre of the condensate, or is weak, then 17 may 
be too small to distort the energy barrier sufficiently, and 
only incomplete cycling occurs even when loj — lOc- This 
feature of the nonlinear system, which agrees with our 
r^ = 3 simulations of the full GP equation, is in contrast 
to the ordinary Rabi case, where complete cycling occurs 
on resonance for any nonzero coupling field. 

The validity of the two-state model breaks down when 
u)f exceeds Wc, because then higher energy vortex eigen 
states are energetically permitted and mixed into the 
state of the system, as seen, for example, in Fig. || and 
the upper dashed curve of Fig. 2. 

The Rabi model provides some insight into the issue 
of the stability of a central vortex state [g 10|. We have 
tested this stability in the T = limit by simulation of 
the GP equation, taking the first excited I — 1 vortex 
state and inserting and withdrawing a narrow stirrer at 
a fixed location in the laboratory frame. We find that al- 
though the condensate then wobbles vigorously, the vor- 
tex undergoes only a very stable small-amplitude oscilla- 
tion about the trap centre |lq| . We can interpret this as 
Rabi cycling of very large effective detuning (i.e. ujf =0), 
and consequently very small transfer probability out of 
the initial vortex state. 

The regular cyclic single-vortex behaviour we have 
found also suggests an experimental technique for prepar- 
ing a condenstate in a central vortex state. By stirring a 
condensate for a half-cycle, a vortex will be drawn into a 
nearly central position. 

In conclusion, we have given a simple, quantitative 
analysis of the single vortex regime of a stirred conden- 
sate. Our two state model captures the essential coherent 
dynamics, and accurately predicts the major features of 
this regime, but also provides a qualitative understanding 
in terms of the concepts of the well known Rabi model. 
Our result for the critical angular frequency can be qual- 
itatively related to that for a rotating cylinder of He II. 
However, in our case the condensate is inhomogeneous, 



and the trapping potential plays a central role, giving rise 
to well separated condensate eigenstates, of which only 
the lowest two become significantly involved. It is worth 
remarking that the speed of sound in the vicinity of the 
perturber has no relevance to the generation of vortices, 
in the scenarios we consider here. The model is also easily 
generalizable to an arbitrarily shaped stirring potential, 
including a rotating anisotropic potential. Our numerical 
calculations have been carried out in two spatial dimen- 
sions, but can be expected to apply to "pancake" con- 
densates, where the dynamics in the axial direction are 
frozen out by very tight axial confinement. Qualitative 
features of our results may have even greater general- 
ity, since the two-state model has no direct dependence 
on dimensionality, and will apply in three dimensions if 
the system symmetry confines the stirrer to couple the 
ground state primarily to a single vortex state. 
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